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Abstract 

A new class of exclusion type processes acting in continuum with synchronous updat- 
ing is introduced and studied. Ergodic averages of particle velocities are obtained and 
their connections to other statistical quantities, in particular to the particle density (the 
so called Fundamental Diagram) is analyzed rigorously. The main technical tool is a 
"dynamical" coupling applied in a nonstandard fashion: we do not prove the existence of 
the successful coupling (which even might not hold) but instead use its presence/absence 
as an important diagnostic tool. Despite that this approach cannot be applied to lattice 
systems directly, it allows to obtain new results for the lattice systems embedding them 
to the systems in continuum. Applications to the traffic flows modelling are discussed as 
well. 

1 Introduction 

The classical simple exclusion process is a Markov process that describes nearest-neighbor ran- 
dom walks of a collection of particles on the one-dimensional infinite 1 integer lattice. Particles 
interact through the hard core exclusion rule, which means that at most one particle is allowed 
at each site. This seemingly very particular process introduced first in 1970 by Frank Spitzer 
[20] appears naturally in a very broad list of scientific fields starting from various models of 
traffic flows [16, 13, 10, 5, 6], molecular motors and protein synthesis in biology(see e.g. [21]), 
surface growth or percolation processes in physics (see [18, 8] for a review), and up to the 
analysis of Young diagrams in Representation Theory [9]. 

Qualitatively from the point of view of the order of particle interactions there are two 
principally different types of exclusion processes: with synchronous and asynchronous updating 
rules. In the latter case at each moment of time a.s. at most one particle may move and hence 
only a single interaction may take place. This is the main model considered in the mathematical 
literature (see e.g. [14, 22, 20, 17] for a general account and [1, 11, 19] for recent results), and 
indeed, the assumption about the asynchronous updating is quite natural in the continuous time 
setting. The synchronous updating means that all particles are trying to move simultaneously 
and hence an arbitrary large (and even infinite) number of interactions may occur at the same 
time. This makes the analysis of the synchronous updating case much more difficult, but this 
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Figure 1: TASEP in continuum. 



is what happens in the discrete time case. 2 This case is much less studied, but still there are a 
few results describing ergodic properties of such processes [3, 5, 6, 7, 10, 13, 16]. 

Our aim is to introduce and study the synchronous updating version of the exclusion process 
in continuum. Note that recently some other interacting particle processes were generalized 
from lattice to continuum case (see e.g. [18, 12]). 

A configuration x := {xi}i £ z is a bi-infinite sequence of real numbers Xi G K. interpreted as 
centers of particles represented by balls of radius r > (see Fig. 1) and ordered with respect 
to their positions (i.e. . . . < X-\ < xq < x\ < . . .). To emphasize the dependence on the radius 
r > we shall use the notation x(r). We say that a configuration x(r) is admissible if 

Xi(r) + r < x i+ i(r) — r Mi G Z 

(the corresponding balls may only touch each other) and denote by X{r) the space of admissible 
configurations. 

The dynamics will be defined as follows. We assume given a collection of (possibly random) 
values {vf}ij, where i,t G Z and t > 0; conditions on this collection will be given shortly. 
For a trivial configuration consisting of a single particle located at time t > at x l G ffi (i.e. 
x l = {xq}) the dynamics is defined as 

Xq := x + v , 

and thus v l Q may be interpreted as a local velocity at time t, i.e. this is simply a random walk on 
M. To generalize this trivial setting for an infinite configuration x(r) G X(r) we again interpret 
a (be-infinite on i G Z) sequence {v\}i )t as local velocities for particles in x*(r) performing 
random walks conditioned to the order preservation and the hard core exclusion rule. 

To simplify presentation we restrict ourselves here to the case of nonnegative local velocities 
postponing the discussion of the general case when the local velocities take both positive and 
negative signs to Section 6. The point is that the formulations in the latter case are becoming 
much more involved, but the results and arguments work with only very slight changes. 

Since only nonnegative local velocities are considered the hard core exclusion rule means 
that the admissibility condition breaks down for the i-th particle at time t G Z + if and only if 
the inequality 

x\{r) +vl + r < 4n(r) -r 

does not hold. If this happens we say that there is a conflict between the particles i and i + 
and to resolve it one applies a normalizing construction 

v \->M{vlx\r)). 

After the normalization the positions of particles are calculated according to the rule 

x l+\r) := x\{r) + M{v\,x\r)) Vi. 
2 if one do not consider some "artificial" updating rules like a sequential or random updating. 
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Figure 2: Positions of particles at time t,t + 1 in cases of weak (a-c) and strong (a'-c') normal- 
izations. Local particle velocities are shown by vectors. The cases (c,c') correspond to negative 
velocities and will be discussed in Section 6. 



In what follows we always assume 3 that Mi,t N(v *, x*(r)) G [0,vj] (to simplify notation by 
the segment [a, b] we mean [min(a, b), max(a, b)]) and consider only nonanticipating normaliza- 
tions 4 satisfying the condition that in the case of the conflict of the i-ih particle with the j-th 
one 5 at time t the position of the i-th particle at the next moment of time :r* +1 (r) G [x*(r), x*-(r)]. 

The normalization may be done in a number of ways and we restrict ourselves to two extreme 
constructions. The first of them we call strong normalization (notation Af s (-, •)) an d according 
to the name we reject (nullify) the velocity leading to the conflict. The second construction we 
call weak normalization (notation J\f w (-, •)) and in this case we modify the conflicting velocity 
to allow the particle to move as far as possible. In terms of gaps 

Aitfir)) = A- := x\ +1 {r) - x\{r) - 2r 

between particles in the configuration x l the normalization procedures are written as follows: 



^^))=={ otherwise ' N w {vl,x\r)) := |^ ^-.^ 



Fig. 2 demonstrates possible positions of particles at two consecutive moments of time t and 
t+l for the cases of weak (a-c) and strong (a'-c') normalizations. Despite appearances these two 
normalization procedures lead to a very different limit behavior of the corresponding particle 
systems. The simplest example (existing even in the continuous time case) is the situation when 
vj = v Vi, t and the gaps between particles in x are smaller than v. Then under the strong 
normalization no motion is allowed, while the weak normalization leads to the well defined 
motion - the exchange of gaps between particles. Other normalization procedures together 
with more general assumptions about the dynamics will be discussed in Section 7. 

Observe that any two particle configurations x(r), x(r) having the same sequence of gaps 
A := {Aj} may be transformed to each other by a one-to-one map 

Xi(r) = <f(xi(r)) := Xi(r) — 2i(r — r) Vi G Z. (1.1) 

Since the normalization procedures that we consider depend only on the gaps between particles 
it is enough to study the case r = 0. On the other hand, if r = 1/2, x®(r) G Z Mi G Z and 



3 This formulation allows to consider velocities of both signs which we shall do in Section 6 and simply means 
that the normalized velocity has the same direction as the original one and cannot exceed it on modulus. 
4 In Section 7 wc shall show that the violation of this condition makes the system to be not well posed. 
5 For nonnegative velocities j = i + 1, but in general j E {i — l,i + l}. 
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v\ G Z Mi G Z, i > then x*(r) 6 Z Vi 6 Z, t > which means that we get a lattice particle 
system. Thus our results lead to a completely new approach to the analysis of lattice systems 
as well. Note however that in the case r = an arbitrary number of particles may share the 
same spatial position which is prohibited in the lattice case. 

Due to the observation above we shall study in detail only the case r = since the corre- 
sponding results for any r > are readily available through the transformation (1.1), see e.g. 
specific calculations for densities and velocities in Lemmas 2.1, 2.4 and Corollaries 4.5, 5.3. 

To simplify notation we shall use the convention x(r) = x°(r), x = x°(0) and similarly 
X = X(0). 

Of course, without some specific assumptions on the structure of local velocities {vj}i it no 
interesting results are possible. We assume that v\ G [0,t>] Mi G Z, t G Z := Z + U {0} and 
one of the following seemingly opposite assumptions holds: 



t-i 




(b) {vj} are i.i.d. (both in i and t) random variables. 



Note that the intersection between the sets of local velocities satisfying the assumptions (a) 
and (b) contains an important case of pure deterministic velocities: v\ = v Mi G Z, t G Z . As 
we shall show properties of systems with local velocities satisfying to the assumption (a) are 
close to the pure deterministic setting. Therefore we refer to the setting (a) as deterministic^ 
and to the setting (b) as random. 

It is of interest that in the seemingly simplest purely deterministic setting v\ = v Mi £ 
Z, t G Z the behavior of the corresponding deterministic dynamical system describing the 
dynamics of particle configurations is far from being trivial. In Section 4.3 we prove that this 
system is chaotic in the sense that its topological entropy is positive (and even infinite). 

To emphasize that under dynamics no creation or annihilation of particles may take place 
this sort of systems is called diffusive driven systems (DDS) instead of a more general object - 
interacting particle systems (IPS). 

The main technical tool in our analysis is a (somewhat unusual) "dynamical" coupling con- 
struction. Despite that various couplings are widely used in the analysis of IPS, applications 
of our approach is very different from conventional. In particular, we do not prove the ex- 
istence of the so called successful coupling (which even might not hold) but instead use its 
presence/absence as an important diagnostic tool. Remark also that typically one uses the 
coupling argument to prove the uniqueness of the invariant measure and to derive later other 
results from this fact. In our case there might be a very large number of ergodic invariant mea- 
sures or no invariant measures at all (recall the trivial example of a single particle performing a 
skewed random walk). The latter example indicates that there is another important statistical 
quantity - average particles velocity that can be computed at least in this case. (See e.g. [2] for 
a discussion of the average velocity in the context of Queueing Networks.) The dynamical cou- 
pling will be used directly to find connections between the average particle velocities and other 
statistical features of the systems under consideration, in particular with the corresponding 
particle densities. 

It is worth note that all approaches used to study lattice versions of DDS are heavily based 
on the combinatorial structure of particle configurations. This structure has no counterparts in 
the continuum setting under consideration. In particular the particle - vacancy symmetry is no 
longer applicable in our case. This explains the need to develop a fundamentally new techniques 

6 In this case Vq might be a trajectory of a deterministic chaotic map / : [0, 1] — ► [0, 1], e.g. v^ +1 := vp^l/v), 
as well as a realization of a true random Markov chain) . 
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for the analysis of DDS in continuum. Despite this new techniques cannot be applied directly 
in the lattice case, the embedding of lattice systems to the continuum setting allows to obtain 
(indirectly) new results for the lattice systems as well. 

The paper is organized as follows. In Section 2 we introduce main statistical quantities under 
study: particle densities, average velocities, etc. and derive their basic properties. Section 3 is 
dedicated to the main technical tool - dynamical coupling. In Section 4 we apply this coupling 
in the weak normalization setting to prove the uniqueness of the average velocity (Theorem 4.1) 
and to derive the complete Fundamental Diagram for the deterministic case (Theorem 4.2). We 
calculate also the topological entropy of this process (Theorem 4.3). The strong normalization 
case is considered in Section 5 (Theorem 5.1), while a more general setting with local velocities 
of both signs is studied in Section 6 (Theorem 6.1). Finally, in Section 7 we discuss some 
generalizations of our results and applications to certain specific traffic models. 

Acknowledgements. The author is grateful to B. Gurevich, S. Pirogov and an anonymous 
referee for a number of valuable remarks. This research has been partially supported by Russian 
Foundation for Fundamental Research, Program ONIT, and French Ministry of Education 
grants. 



2 Basic properties of DDS 

Here we shall study questions related to densities and velocities of DDS. To simplify notation we 
use the convention that the normalization N G {Af s ,Af w } and specify it only if this is necessary. 

By the density p(x, I) of a configuration x G X in a bounded segment I — [a, b] G M we 
mean the number of particles from x whose centers X{ belong to I divided by the Lebesgue 
measure |/| > of the segment /. If for any sequence of nested bounded segments {I n } with 



|/„| ""l^ qq the limit 



p(x) := lim p(x, I„ 



exists and does not depend on {/„} we call it the density 7 of the configuration x G X . Otherwise 
one considers upper and lower particle densities p±{x) corresponding to upper and lower limits. 

The correspondence between particle densities for configurations with r = and r > is 
given by the following statement. 

Lemma 2.1 Let configurations x(r) G X(r), r > and x G X have the same sequence of gaps 
{A,}. Thenp ± (x(r)) = T J^. 

Proof. Due to the one-to-one correspondence (1.1) between the configurations x(r) and x, 
for each segment J C 1R 1 which contains p(x,I) ■ \I\ particles from the configuration x, one 
constructs the segment I(r) containing the same particles from the configuration x(r). The 
length of this segment is equal to \I(r)\ = \I\ + 2r • p(x, I) ■ \I\. Therefore 

p(x(r),I(r)) - 



\I\ +2rp(x,I) ■ \I\ l + 2rp(x,I)' 
Passing to the limit as |/| — > oo one gets the result. □ 



Remark 2.2 If 3p(x) < oo then \x n — x m \/\n — m\ — > p{%)- 

7 In Section 7.2 we shall show that this definition may be significantly weaken in the case when all particles 
move in the same direction. 
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Lemma 2.3 The upper/lower densities p±{x t ) are preserved by dynamics, i.e. p±(x t ) = 

P± {x t+1 ) vtez . 

Proof. For a given segment / £ R the number of particles from the configuration x* G X 
which can leave it during the next time step cannot exceed 1 and the number of particles which 
can enter this segment also cannot exceed 1. Thus the total change of the number of particles 
in / cannot exceed 1, because if a particle leaves the segment through one of its ends no other 
particle can enter through this end. Therefore 

\p{x\l)-p{x t+ \l)\.\I\<\ 

which implies the claim. □ 

By the (average) velocity of the i-th particle in the configuration x G X at time t > we 
mean 

V(x,i,t) := j ^A/"(^,:r s ) = (xl-x^/t. 

If the limit 

V(x,i) : = lim V(x, i,t) 

exists we call it the (average) velocity of the i-th particle. Otherwise one considers upper and 
lower particle velocities V±(x,i). 

The correspondence between average particle velocities for configurations with r = and 
r > is even simpler than for densities. 

Lemma 2.4 Let configurations x(r) G X(r), r > and x G X have the same sequence of gaps 
{Aj}. Then Wi,t V(x(r),i,t) = V(x,i,t) for a given collection of local velocities {f*}i,t. 

Proof. Observe that the motion of particles depends only on the local velocities and the 
sequence of gaps. Thus at any time t > the sequence of gaps being changing in time is still 
the same for both configurations x(r) and x. Therefore 

^K,x*(r))=jV(t;?,x*) Vi,f 
which yields the claim. □ 



Lemma 2.5 Let x G X then \V(x,j, t) — V(x, i,t) \ * z ^ > a.s. Vi,j G Z. 

Proof. It is enough to prove this result for j = Consider the difference between (average) 
velocities of consecutive particles 

V(x,i + l,t)-V(x,i,t) = ^±l_3±i_^_^ 

X l+l ~ X i _ ~ x i 

= A*//- A°/t. 

The last term vanishes as t — > oo and it is enough to show that the same happens with Aj/t. 
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Consider first the deterministic setting (i.e. vj = Vq) and show that 8 V«, t 

t fmax(v, A°) if AT = AC 

|max(2t>, A?) ii M — M s " 1 J 

Obviously this is true for t = 0. Assume that this inequality holds up to time t G Z and 
consider the moment £ + 1. There might be two two possibilities: 

(a) A- > Vq. Then M(vl : x f ) = Vq and 

A< +1 = A* - M(vj, x«) + M^x*) < A* - v$ + v$ = A* < max^, A°) 
by the assumption. 

(b) A^ < ^. Then Af^v^x 1 ) = A* and A/" s (t>£,£*) = 0. Therefore 

A< +1 = A< - A< + M(vj +1 , x l )<v< max^, A?) if M = M w , 
AI+ 1 =A l t -0 + M(vj +1 , x 1 ) <2v if M = M s . 

Thus in the deterministic setting the gaps are uniformly bounded in time and hence A\/t t -—^> 

0. 

Analysis of the random setting is much more involved since the gaps between particles in 
principle may grow with time and become arbitrary large but this may happen only very slowly. 
To estimate from above the value of the i-th gap A* we drop from the consideration all particles 
except the i-th and (i + l)-th (preserving for all i G Z the velocities {vj, vj +1 } t ) and denote 
the resulting configuration by x l := {x\, x\ +1 } and the gap between this pair of particles by A*. 
We have 

-^Al-Nivlx^+Nivl^l 



A* +1 := A* - M(vj, x<) + M(vj +1 , **) = A* - M (vj, *«) + 



i+i- 



The comparison between A* and A* will be done by induction separately for the weak and 
strong normalizations. 

First let us prove that A* > A* if M = M w . At time t = obviously A° = A°. Assume 
that A* > A* for some t G Z + . Clearly, 

0<^ +1 ,x*)<^ +1 . 

For vj there might be two possibilities: 

(a) vj < A*. Then Mivj^x 1 ) = A/"(u*, x*) = v\ and hence 

A< +1 = A\-v\ + vj +1 > A* - v\ + vj +1 = A t+1 . 

(b) v\ > A*. Then M w (vj, x f ) = A*, M w (vj, > A* and hence 

A* +1 = A • - M w (vj, **) + > = 

If A/" = M s a weaker estimate A* + v > A* takes place. Considering again the same 
possibilities we see that the cases t — and (a) hold without any changes, but the case (b) 
should be rewritten. 



i li Vq takes both positive and negative values then A* < max(4w, A?). 
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(b') v\ > A*. Then A/" s (^,x*) = 0, N a {v\,&) 
Afsiylx*). Thus 

A* +1 =Aj-Af s (vj,x t )+v t l+1 

Consider now the behavior of A* as a function of time t. If A* > v we get vj < A* 
and hence A/^uf, x*) = which implies that outside of the region [0,v] the sequence A* 
behave as a spatially homogeneous reflected at random walk with i.i.d. symmetric increments 
vj +1 —v\. Thus the mathematical expectation -E(A-) cannot exceed 9 2v and hence by Chebyshev 
inequality the probability 

P{K\/t >e)<\ E(Aj/t) <|^0, 
which finishes the proof. □ 



A* if 2 < A* ' and hence ^^'^ - 



- (Ar s (^^)-AT s (^5 t ))>A* +1 -t;. 



Corollary 2.6 TTie upper and lower particle velocities V±(x,i) do not depend on i (but might 
be random). 

3 Coupling 

Recall that a coupling of two Markov chains x* and y* acting on the space X is an arrangement 
of a pair of processes on a common probability space to facilitate their direct comparison, 
namely this is a pairs process defined on the direct product space X x X satisfying the 

assumptions 

P{{x\ y t )eAxX)= P(x t G A) and P((x\ y l ) G X x A) = P(y l G A) 

for any measurable subset iCI, i.e. the projections behave as the individual processes. 

Let x f , x l be two copies of Markov chains, describing the DDS which we consider throughout 
the paper. Typically in continuous time interacting lattice particle systems one uses (see e.g. 
[14]) an equal coupling (pairing) when particles sharing the same sites in the copies are 
considered to be paired and all choices of their velocities are assumed to be identical. This 
sort of coupling works rather well for continuous time systems when only a single particle 
may move at a given moment of time. In the discrete time case the situation is much more 
complicated since an arbitrary number of particles may move simultaneously and thus it is 
possible that the particles of the processes x l \x* pass each other and never share the same 
positions. In fact, this difficulty is not really crucial and can be cured under some simple 
technical assumptions. A more important obstacle is that if a pair is created and only one 
of its members is blocked by an unpaired particle, then due to the simultaneous motion of 
the blocking unpaired particle and the non-blocked particle belonging to the pair the following 
situation may happen: *° — >■ ° °. Thus the old pair will be destroyed but no new pair will 
be created under the equal pairing construction. Here and in the sequel we use a diagrammatic 
representation for coupled configurations, where paired particles are denoted by black circles 

9 4w if local velocities take both positive and negative values. 
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and unpaired ones by open circles, and use the upper line of the diagram for the x-particles 
(i.e. particles from the x-process) and the lower line for the x-particles. 

To deal with this obstacle we introduce a dynamical 10 coupling, a very preliminary version 
of which was described in [7] for the lattice case and was inspired by the idea proposed by 
L. Gray for the simplest discrete time lattice TASEP (unpublished). It is worth mention also 
the coupling proposed for the lattice continuous time case by O. Angel (see [1, 11]). As we 
shall show an important advantage of the dynamical coupling with respect to the Angel's 
construction is that the former guarantees that the distances between mutually paired particles 
are uniformly bounded. 11 

By the dynamical coupling of the processes x f , x l we mean a gradual pairing of close enough 
particles belonging to the opposite processes satisfying the following assumptions: 

(Al) At t — all particles are assumed to be unpaired. Velocities of mutually paired particles 
are identical. 

(A2) Once being created a pair of particles remains present 12 for any moment of time in the 
future, however at different moments of time the roles of the pair's members may be 
played by different particles. 

(A3) A particle overtaking during one time step of the dynamics some unpaired particles from 
the opposite process becomes paired with one of them. 

According to (Al)-(A3) particles from the same pair move synchronously until either the 
admissibility condition breaks down for only one of the particles (which means that its move- 
ment is blocked by another particle) or one of the members of the pair is swapped with an 
unpaired particle from the same process (see Fig. 3 for the case of the weak normalization). It 
is convenient to think about the coupled process as a "gas" of single (unpaired) particles and 
"dumbbells" (pairs). A previously paired particle may inherit the role of the unpaired one from 
one of its neighbors. In order to keep track of positions of unpaired particles we shall refer to 
them as x- and x- defects depending on the process they belong. 

There are a number of ways to realize the dynamical coupling (in particular, using only the 
idea of the particle's overtaking). To demonstrate the flexibility of our approach we describe a 
different construction. Note that in the sequel we shall use only the properties (Al)-(A3) and 
the proofs will not depend on other details of the coupling. 

By the x-triple ( ,° * or * °. ) in the coupled process {x t ,x t ') we mean two mutually paired 
particles and a x-defect located in the segment between them, whose index differs by one from 
the index of the paired x-particle. The f-triple ( * G . or . G *) is defined similarly. 

Two pairs of particles are said to cross each other if straight lines connecting positions of 
particles belonging to the same pair intersect, e.g. % *. , where particles belonging to the same 
pair are marked similarly. 

A rr-defect at x\ together with the closest 13 i-defect at x * ( Q ° or ° Q ) are said to be a d-pair 
if \x\ — Xj\ < v, this pair of defects does not cross with any mutually paired particles, and the 
open segment does not contain any other defects. We say that a d-pair is smaller 

10 The word "dynamical" is meant to emphasize that the mutual arrangement of particles in pairs may change 
with time under dynamics in distinction to the conventual equal coupling (where the particles have coinciding 
positions). 

11 In the Angel's construction the distances may grow to infinity. 

12 Starting from the moment when a pair is created we consider it as an entity independently on the possible 
change of particles forming it. 

13 If there are several closest x -defects one chooses the defect with the smallest index. 
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Figure 3: Pairing of particles. Black circles corresponds to paired particles and open circles to 
defects. The paired particles are connected by straight lines. At time t the particles i and j 
are paired, while at time t + 1 the x-particle i becomes unpaired and the x-particle j becomes 
paired with the re-particle i + 1. The unpaired initially particles i + 2 and j + 1 become paired 
at time t+1. 



than a d-pair (n,m) if \i\ < \n\, or if % < n in case \i\ = \n\. Observe that i = n but j ^ m 
cannot happen in distinction to i ^ n but j = m. 

Note that in the collection °* * the first two re-particles together with the first x-particle 
form a x-triple despite the presence of an additional paired particle in the segment between 
them. On the other hand, the collection * # ° does not contain neither triples nor d-pairs. 

A pair of configurations (x*,x*) representing the coupled process at time t is said to be 
proper if it does not contain x- or x-triples, d-pairs, and crossing mutually paired particles. 

The fact that at time t the pair of configurations (x*,x*) were proper does not imply that 
it remains proper under dynamics at time (t + 1). In particular, triples of both types and 
d-pairs may be created, e.g. , * G — G * or ° G0 — ► ° , however due to the particle order 
preservation crossing mutually paired particles cannot appear. 

Lemma 3.1 Let a pair of configurations (x*,x*) have no crossing mutually paired particles. 
Then among triples of the same kind there are no common elements. 

Proof. Direct inspection. As an illustration let us check the claim about x-triples. Assume 
that two x-triples have a common x-defect (mutually paired particles cannot be common by 
definition). Then this implies that the mutually paired particles in these triples either cross 
each other * + ° * or the index of one of the paired x-particles differs from the index of the 
common defect by more than one * ' ° + , . The latter contradicts to the definition of the 
x-triple, why the former contradicts to the assumption about the absence of crossing mutually 
paired particles. In the diagrams above paired particles from the 2nd triple are marked by stars 
to distinguish them from the 1st triple. □ 



Therefore all triples of the same kind may be resolved simultaneously. This will be done as 
follows. A x- or x-triple is transformed such that the former defect is becoming paired to the 
particle from another process, while another previously paired particle is becoming unpaired: 

o • • o 

The case of a d-pair is even simpler, namely the defects "annihilate" forming mutually 
paired particles: Q ° — >■ .* . In all cases the positions of particles are preserved but their 
"roles" are changing. 

Finally the coupling procedure consists of the following steps: 

(1) Each x-triple is recursively resolved: # ° * — > .* ° . 
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(2) Each x-triple is recursively resolved: * Q . — > % . 

(3) The smallest 14 d-pair is recursively resolved: Q ° — > / . 

Lemma 3.2 The coupling procedure described above is well defined, leads to the Markovian 
coupling, and satisfies the assumptions (A1)-(A3). 



Proof. Let us check that this procedure is well defined. By Lemma 3.1 if a particle belongs 
to a certain triple then it cannot belong to any other triple. On the other hand, segments 
belonging to paired particles may overlap and resolving a x- or f-triple one may create a new- 
one of the same kind: 

This explains the necessity of the recursion during the first two steps of the procedure. Note 
that resolving a x-triple one cannot create a new rr-triple and vice versa (defects do not move 
from one process to another). 

Elements of the smallest d-pair might belong to some other d-pairs. Therefore resolving 
it we might change the d-order of the remaining d-pairs. To take this into account we are 
recalculating the d-order after each recursion procedure. 

Consider now the motion of a given defect under the recursions in the coupling procedure. 
Observe that the defect may move arbitrary far in any direction from its initial position due to 
these recursions: 

o • ••••••• o" 

Nevertheless a defect cannot change its direction of movement. Assume from the contrary that 
a x-defect during two consequent steps of the recursion moved first to the right ( # ° * — ► / °) 
and then to the left (* °. — > ° *.). This can happen only if after the first step of the recursion 
the defect became a member of a new a>triple of type # ° *. Then the only candidate for the 
role of the paired x-particle in this x-triple is the paired x-particle which played the role of this 
defect on the previous recursion step. We came to the contradiction, because a particle may 
belong to only one pair. 

Thus the recursion is finite in the sense that each defect in a bounded spatial segment in 
finite time either will stop moving or will leave this segment and never return back. Note 
however that in general one cannot divide a configuration into finite pieces and deal with them 
separately since a defect may move from one piece to another. 

After the application of the first two steps all x- or i-triples will be eliminated and only 
d-pairs may be present. Observe now that when one resolves a d-pair neither triples nor 
new defects are created. However since various d-pairs may intersect they should be resolved 
separately during the last step. Additionally neither of above procedures may create crossing 
pairs of mutually paired particles (since members of different triples of the same type do not 
intersect and c- and d-pairs cannot cross each other). 

Let the pair of configurations (x*" 1 , a;*" 1 ) be proper. Then according to arguments above 
after one time step of the dynamics the application of the coupling procedure, is well defined 
and the pair of configurations at time t is proper as well. 

By the construction the one-time step transition probabilities for both processes x f and x l 
remain unchanged and the one-time step transition probabilities for the pairs process are well 
defined. Therefore this construction defines a Markovian coupling between two copies of the 
Markov chain describing our DDS. 

14 The ordering of d-pairs is updated after each recursion procedure. 
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The property (Al) holds by the construction. A pair breaks down only if one of its members 
is replaced by an unpaired particle, and hence the pair as a whole survives. This proves (A2). 
The property (A3) follows from the fact that under the one time step of the dynamics of a 
proper pair of configurations all objects under consideration: x- and f-triples, and d-pairs may 
be created only during the particles overtaking. □ 

Denote by p u (x,I) the density of the x-defects belonging to a finite segment J, and by 
p u (x) := p u (x, R) the upper limit of p u (x, I n ) taken over all possible collections of nested finite 
segments I n whose lengths go to infinity. 

We say that a coupling of two Markov particle processes is nearly successful if the 

upper density of the x-defects p u (x) vanishes with time a.s. This definition differs significantly 
from the conventional definition of the successful coupling (see e.e. [14]), which basically means 
that the coupled processes converge to each other in finite time. 

In the random setting under some regularity assumptions the dynamical coupling turns out 
to be nearly successful (the proof of this result goes out of the scope of the present paper 
and will be published elsewhere), however in general especially in the deterministic setting this 
property needs not hold. 

Applying the notion of the nearly successful coupling to the exclusion process under study 
we get the following conditional result. 

Lemma 3.3 Let x, x G X with p(x) = p(x), and let there exist a nearly successful coupling 
(x*,f*) such that distances between the pair members are uniformly bounded from above by 
7(t) = o(t). Then 

\V(x,0,t) -V(x,0,t)\ *-=^0. 

Proof. Consider an integer valued function n t which is equal to the index of the x-particle 

paired at time t > with the 0-th x-particle. If the 0-th x-particle is not paired at time t we 

frii-i if t > 
setn,:=| if t = 

To estimate the growth rate of \n t \ at large t observe that n t changes its value only at those 
moments of time when the 0-th x-particle meets a rr-defect. By the assumption about the 
nearly successful coupling at time t ^> 1 the average distance between the defects at time t is 
of order 1/ p u (x l ) while the amount of time needed for two particles separated by the distance 
L to meet cannot be smaller than L/{2v). Therefore the frequency of interactions of the 0-th 
re-particle with f-defects may be estimated from above by the quantity of order p u (x f ) 0, 
which implies n t /t t —^+ 0. 

Now we are ready to prove the main claim. 

\V(x,0,t)-V(x,0,t)\ = \(xl-x° ) - (4 -x° )\/t 

< \x — x \/t + \x — x \/t 

< |4 - 4j/f + M \± t nt - 41/hi + K -x° \/t. 

The 1st addend can be estimated from above by i{t)/t 0. The 2nd addend is a product of 
two terms \n t \/t and \x l nt — x\\/\n t \. As we have shown, the 1st of them vanishes with time. If 
\n t \ is uniformly bounded, then the 2nd term is obviously uniformly bounded on t. Otherwise, 
for large \n t \ by Remark 2.2 and the density preservation the 2nd term is of order p(x), which 
proves its uniform boundedness as well. Thus the 2nd addend goes to as t — > oo. Noting 
finally that the last addend also vanishes with time we are getting the result. □ 
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4 Weak normalization 



Consider the coupled process under the weak normalization and set Wfj := x\ — x*. 

Lemma 4.1 The supremum of \Wfj\ taken over all mutually paired particles is uniformly 
bounded by v for any t G Z . 

Proof. We start at time t = when there are no pairs and wait until the first of them appears. 
At that moment the distance between the members in a pair cannot exceed v. Starting from 
that moment the distances may grow and some new pairs may be created. Contrary to our 
claim assume that there is the first moment of time t at which there is a pair of particles located 
at x\,x^ for which \x\ — > v and it is the largest distance between the paired particles at 
that moment of time (or one of the largest) and such that — x t ~ 1 \ < v. According to 
the definition of the pairing process there are no unpaired particles between the particles from 
the same pair. Therefore in order to enlarge the distance between the particles one of them 
should be blocked by a particle from another pair, which contradicts to the assumption about 
the maximality of the distance. □ 



Lemma 4.2 Let p(x) = p{x) and let in the coupled process 3 a (random) moment of 

time tij < oo such that x\ > f* for each t > Uj. Then the coupling is nearly successful. 

Proof. By the assumption each x-particle will overtake eventually each x-particle located 
originally to the right from its own position and thus will form a pair with it or with one of its 
neighbors (if they are so close that were overtaken simultaneously). Thus the creation of pairs 
is unavoidable. To show that the upper density of defects cannot remain positive, consider 
how the defects move under our assumptions. Assume that at time t > the i-th x-particle 
is paired with the j-th x -particle. Then by Lemma 4.1 in order to overtake at time s > t the 
j-th re-particle significantly (by a distance larger than v ) the i-th rc-particle necessarily needs 
to break the pairing with the j-th x-particle. Thus by the property (A3) of the dynamical 
coupling either a x-defect overtakes the j-th i-particle: °. * — > ,° * — >■ * °, or the 
i-th x-particle overtakes a x -defect: . * — > , ' — >■ .*. (Otherwise this pair will not be 
broken.) Therefore during this process the x-defects move to the right while the x-defects move 
to the left. Hence they inevitably meet each other and "annihilate". The assumption about 
the equality of particle densities implies the result. □ 



4.1 Uniqueness of the average velocity 

As we shall see under our assumptions even in the weak normalization case the nearly successful 
coupling needs not hold (e.g. in the deterministic setting). Therefore one cannot apply directly 
Lemma 3.3 in this case. Nevertheless we shall show that the absence of coupling is not a serious 
obstacle and it can be used as a diagnostic tool. 

Theorem 4.1 In the weak normalization case the set of limit points as t — > oo of the sequence 
{V(x, t)} t £z depends only on the density p(x) assuming that the latter is well defined. 

Proof. Consider a general DDS under the weak normalization. Let x, x G X p := {z G 
X : p(z) = p} be two admissible configurations of the same particle density. If one assumes 
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that the coupling procedure described in Section 3 leads to the nearly successful coupling of 
particles in these configurations then by Lemma 4.1 the assumptions of Lemma 3.3 are satisfied 
and hence \V(x, 0, t) — V(x, 0, t) \ < -^5 which by Lemma 2.5 implies the claim. In general the 
assumption about the nearly successful coupling may not hold, 15 however as we demonstrate 
below the pairing construction is still applicable. 
Define random variables 

Wj. :=x\-x), i,jeZ, teZ - 

Then 

V(x,i,t)-V(xJ,t) = Wt j /t-W? j /t. 

Since by Lemma 2.5 the differences between average velocities of different particles belonging 
to the same configuration vanish with time it is enough to consider only the case % — j — 0. 
For Woo there might be three possibilities which we study separately: 

(a) lim W^/t = 0. Then \V(x, 0, t) - V(x, 0, t)\ < \W^\/t + \W^\/t *-=^> 0, which by Corol- 
lary 2.6 implies that the sets of limit points of the average velocities coincide. 

(b) lim sup W^/t > 0. Then Mi E Z the i-th particle of the x-process will overtake eventually 

t— *oo 

each particle of the i-process located at time t — to the right from the point x®. 
This together with the assumption of the equality of particle densities allows to apply 
Lemma 4.2 according to which the coupling is nearly successful. On the other hand, by 
Lemma 4.1 the distance between mutually paired particles cannot exceed v. Therefore by 
Lemma 3.3 we have \V(x, 0,t) — V(x, 0,t) \ 0, which contradicts to the assumption 
(b). 

(c) lim sup W^oo/t < 0. Changing the roles of the processes x f , x l one reduces this case to the 
case (b). 

Thus only the case (a) may take place. □ 



4.2 Deterministic setting 

Theorem 4.2 (Fundamental Diagram) In the deterministic setting 

V( X) = Ita igmWl/pX) = {l /p{x) «J™ ^ UV (4.1) 

ifv f = v. 

Proof. Consider a family 

X p := {x E X : Xi := i/p + co, cu E M} 

of uniformly spatially distributed configurations of a given density p > 0. This set is forward 
invariant and 

xl +l -x\ = min(l/p, u*) \fx l E X p ,i E Z, 

15 Consider e.g. the deterministic setting with 1/p > 5v and the configurations Xi :— i/p and Xi := i/p+ 2v. 
Then p(x) — p(x) — p, V(x) — V(x) — v but no pair will be created. 
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i.e. all particles in the configuration get the same normalized local velocity min(l/p, Vq) (de- 
pending in general on time t). By the definition of the deterministic setting the limit 

1 

V(x) := lim - Vmin(l/p,^) 
t— +00 t ' 

is well defined. On the other hand, by Theorem 4.1 all configurations of the same density have 
the same average velocity, which implies the result. □ 



Remark 4.3 This result looks very similar to the one known for the deterministic version 
of the lattice TASEP (see [16, 5]), however the latter case is characterized by the following 
feature: if the density is large enough particles inevitably form dense clusters without vacancies 
inside (static traffic jams). The proof above shows that the "typical" behavior of high density 
configurations in continuum is different: they do form particle clusters, but these clusters are 
not staying at rest but are moving at a constant velocity as an "echelon" . It is of interest that 
in order to imitate such behavior a number of complicated lattice models were developed. 

Remark 4.4 The construction used in the proof is especially striking in that the same family 
of uniformly spatially distributed configurations allows to study the limit dynamics in the 
deterministic setting for all configurations having densities. Note that this argument cannot be 
applied directly in the lattice version of DDS. Nevertheless since the "lattice configurations" 
are included in DDS under consideration the result holds as well, which implies completely new 
results for lattice TASEPs with long jumps. 

Corollary 4.5 Let x(r) G X(r), r > and p(x(r)) be well defined and let \ti,t vj = v. Then 

, (v if p(x) < — ^- 

V(x(r)) = \ \u v+2r ■ 

I l/p{x{r)) — 2r otherwise 



In particular in the lattice setting this reads 



V(x(l/2)) 



v 



if P(*)<4 



l/p(x(l/2)) — 1 otherwise 



+1 



Proof. By (1.1) and Lemma 2.1 for each configuration x(r) one constructs the configuration 
x with the same sequence of gaps and the relation between their densities is written as 

p{x) - P( " (r)) 



1 — 2rp(x(r)) 



Additionally by Lemma 2.4 average velocities related to configurations with the same sequence 
of gaps coincide. Substituting p(x) as a function of p(x(r)) to (4.1) we get the result. □ 



4.3 Entropy 

In this Section we restrict the analysis to the pure deterministic setting (i.e. v\ = v Vi,t). 
Then our DDS is defined by a deterministic map T v : X — > X from the set of admissible 
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configurations into itself. Our aim is to show that this map is chaotic in the sense that its 
topological entropy is infinite. 16 

We refer the reader to [4, 23] for detailed definitions of the topological and metric entropies 
for deterministic dynamical systems and their properties that we use here. To avoid difficulties 
related to the non-compactness of the phase space we define the topological entropy of a map T v 
(notation h top (T v )) as the supremum of metric entropies of this map taken over all probabilistic 
invariant measures (compare to the conventional definition of the topological entropy and its 
properties in [23]). 

For a finite subset of integers / and a collection C := {Ci}i e j of open intervals the subset 
Ci y c ■= {x £ X : Xi G Ci \/i G /} is called a finite cylinder. 17 We endow the space of 
admissible configurations X by the a-algebra B generated by the finite cylinders defining a 
topology in this space. 

We start the analysis with the action of a shift-map in continuum a v : X — > X defined as 

(cr v x)i := Xi + v i G Z,x G X. 
Lemma 4.6 The topological entropy of the shift-map in continuum a v is infinite. 

Proof. The preimage of a finite cylinder under the action of a v is again a finite cylinder. 
Therefore this map is continuous in the topology induced by the cr-algebra B generated by 
finite cylinders. 

The idea of the proof is to construct an invariant subset of X on which the map o~ v is 
isomorphic to the full shift-map in the space of sequences with a countable alphabet. The 
result follows from the observation that the topological entropy of the full shift-map with 
the alphabet consisting of n elements is equal to Inn (see, e.g. [4, 23]). 

Let a := {ai} iG z + with ctj G (0,t>) and let a n := {aj}™ =1 . Consider a sequence of subsets 
X*™) c X consisting of all configurations x G X satisfying the condition Wk G Z X2k G 
vZ, X2k+i G x 2 k + ot n '. Then X^ is ovinvariant and the restriction a v \X^ is isomorphic 
to the full shift-map with the alphabet A n consisting of n elements {a^} of type : = 
{[0,aj), [ai,v)}, i.e. each element is represented by a pair of neighboring intervals. Therefore 
the topological entropy of is equal to Inn "— ? oo. □ 

Another elegant (but technically difficult) way to derive this result was proposed by Boris 
Gurevich. Consider a special flow S* corresponding to the shift-map acting on the sequences 
{Ai(x)} with the roof function equal to the first nonnegative particle coordinate. This shift-map 
has an infinite alphabet, hence its entropy is infinite. The special flow S 1 is isomorphic to the 
1-shift of {x^, while the entropy of the special flow can be calculated by the Abramov-Rohlin 
formula. 

Theorem 4.3 The topological entropy of the pure deterministic exclusion process in continuum 
is infinite. 

Proof. The preimage of a finite cylinder under the action of T v is again a finite cylinder. 
Therefore this map is continuous in the topology induced by the a-algebra B generated by 
finite cylinders. 

16 Normally one says that a map is chaotic if its topological entropy is positive, so infinite value of the entropy 
indicates a very high level of chaoticity. 

17 In general the cylinder C/,c might be empty for nonempty sets /, C. 
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Observe that the subset X := {x e X : A;(x) > t> Vi G Z} of the set of admissible 
configurations is T„-invariant. Therefore h top (T v ) > h top (T v \X ) and for our purposes it is 
enough to show that the latter is infinite. On the other hand, by the definition of the map T v 
we have T V \X = a v \X . 

We still cannot apply the result of Lemma 4.6 directly because in the case under consider- 
ation the gaps between particles are greater or equal to v by the construction. Recall that in 
the proof of Lemma 4.6 the gaps were not greater than v. To this end one sets ctj G (v, 2v) and 
modifies the definition of X^ as follows: 

x 2 k+i e x 2 k + a n VA; e Z, z 2 fc e 3uZ. 

Consider the the alphabet with elements of type := {[0,aj), [aj,3f)}. Then the 3-d 
power of the map T v \Xq is isomorphic to the full shift-map with the alphabet A^ n \ Using 
that 

3h top (T v \X ) = h top ((T v \X f) = h top (aW) = Inn 
we get the result. □ 



5 Strong normalization 

Recall that W-j := x\ - x * for x\ x l e X, t>0. 

Lemma 5.1 There exists a coupled process such that under the strong normalization 

supjj t Wfj = oo, where the supremum is taken over all mutually paired particles. 

Proof. It seems that the argument applied in the weak normalization case should work 
also in the case of the strong normalization. However, a close look shows that in this case a 
"blocked" particle does not move to "touch" the particle conflicting with it (as it would in the 
weak normalization case) but preserves its position instead. Therefore the distance between 
members of the same pair may become larger than the distance between the members of the 
"blocking" pair which cannot happen in the weak normalization case: . *. " — > . ** *. 
Here initially distances between members in pairs do not exceed v. The 1st pair is blocked by 
the 2nd pair and since the i-member of the 1st pair cannot move (while the x-member can) 
the distance between them becomes larger than v. 

To demonstrate that distances between members in pairs may grow to infinity fix some 
< e 1 and consider a pair of configurations x, x such that Xq = Xq = and A 2 k = 
|(i> — e), A 2 k+i = \{v — s), Afc = v — e VA; e Z. After the application of the pairing procedure 
V? the i-particles in both configurations will become paired forever. On the other hand, under 
dynamics x l = x° Vi while the x-particles having gaps greater than v will at constant velocity 
v. Therefore the distances between members in pairs will grow linearly with time. □ 

This result demonstrates and partially explains a significant difference in the behavior of 
DDS under weak and strong normalizations. Still, as we are going to show, at least some 
features of the Fundamental Diagram are preserved. Consider the pure deterministic setting 
(i.e. v\ = v). The inequality (2.1) shows that in this case gaps between particles cannot 
become much larger than their initial values. The following result demonstrates that under 
some mild additional assumptions (which definitely hold for high particle densities) large gaps 
will disappear with time. 
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Figure 4: Fundamental Diagram (dependence of the average velocity V on the particle density 
p) for the pure deterministic setting under the strong normalization. The curvilinear region 
H :— {(p, V) : - p — v <V <- p , V < v} corresponds to the hysteresis phase. 



Lemma 5.2 Let x G X be spatially periodic and we consider only the pure deterministic setting 
(i.e. v\ = v). Assume that Vt 3j > t : A j (x t ) < v. Then Vi 3t { < oo : A^x*) < 2v \/t > t { . 

Proof. Observe that the spatial periodicity and its period is preserved under the pure de- 
terministic dynamics. Thus the situation is equivalent to the consideration of a finite number 
(say N) particles on a ring and to the assumption that for each t G Z + among these particles 
there is a particle with a gap less than v ahead of it. Note that according to the definition of 
the strong normalization J\T s {v\,x t ') = whenever Aj(x*) < v. By (2.1) Aj(x*) < 2v implies 
Ai(x t+1 ) < 2v. Therefore new new long gaps (of size larger or equal to 2v) cannot be created 
and we need to show only that long gaps in the original configuration will cease to exist with 
time. 

By the assumption for any t there exists a short gap (shorter than v) and the corresponding 
particle will not move during the next time step. Thus the index of the short gap decreases by 
one after each time step until it "collides" with one of the long gaps: Aj(x*) > 2v, A i+1 (x l ) < v. 
On the next time step A;(x m ) := Aj(x*) — v. Due the spatial periodicity the amount of time 
between these "collisions" is bounded and after each of them the length of a long gap decreases 
by v. Thus they will disappear in finite time. □ 



Theorem 5.1 Let x E X and p(x) be well defined. Then V(x) = v if p(x) < ^ and otherwise 
for a.e. point (p,V) in the curvilinear region 

H := {(p,V) : max(l/p-v,0) < V < min(l/p,u)} 

(see Fig. 4) there exists a configuration x G X with p(x) = p,V(x) = V, i.e. the region H 
corresponds to the hysteresis. 

Proof. We say that particles numbered from i + 1 to i + k with i G Z, k G Z + belonging 
to an admissible configuration x G X form a cluster of length k if all gaps between them 
are strictly less than v and the gaps to surrounding particles are not smaller than v, i.e. 
A i+ j < v \fj — 1, 2, . . . , k — 1 and Aj, A i+k > v. Positions of particles belonging to the cluster 
are changing with time, and leading particles leave it, while some new particles may join the 
cluster from the other side. Nevertheless the length of a cluster cannot grow with time (and 
new clusters cannot be born in the pure deterministic setting in distinction to the random one) 
since the rate with which the leading particle leaves the cluster (one per unit time) is at least 
not smaller than the rate at which new particles join the cluster from the other side. 
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We start with the analysis of configurations of low density (smaller than ^) and our aim 
is to show that in this case each particle achieves eventually the largest available velocity v. 
Consider the motion of the 0-th particle in a configuration x G X with < p(x) < j- and 
denote by t the first moment of time after which this particle will not join any cluster. If t < oo 
then Afsiyl) = Vt > i and hence V(x, 0, t) v. 

If t — oo then there exists an infinite sequence of clusters of growing length such that the 
0-th particle joins each of them consecutively. Let us show that this assumption contradicts 
to the condition that p(x) < We number the clusters to which the 0-th particle will join 
according to their natural order starting from k = 1 and introduce the following notation: t k - 
the moment of time when the 0-th particle joins the fc-th cluster, n k - the number of particles 
in this cluster, m k - the number of particles in the open segment between xq and the beginning 
of this cluster, and L k - the length of the minimal segment containing the fc-th cluster and the 
point x . Then 

/ / r i\ m k ~\~ k— too , \ 

p{x, {x ,x + L k \) = ► p(x). 

J^k 

All m k particles will join the fc-th cluster during the time t k and at time t k this cluster should 
still exist. Therefore the distance which the 0-th particle covers during this time cannot be 
smaller than L k — m k v — n k v while its velocity cannot exceed v and thus 

t k v > L k - m k v - n k v. 

On the other hand, exactly t k particles will leave the cluster during this time, i.e. m k + n k > t k . 
This gives 

m k + n k > tk_ > L k /v -m k -n k = I _ m k + n k ^ 
Lk ~~ L k ~~ L k v L k 



Therefore 



- < 2 > 2p(x), 

v L k 



which proves the desired claim that t = oo implies p(x) > 

Consider now the case of densities greater than ^. In this case there might be two possi- 
bilities: 

(a) All particles will eventually achieve the largest available velocity v. Then the gaps will 
become not smaller than v and hence they cannot exceed 2v (by the assumption on the density 
region). Obviously this situation may take place only if p{x) G [^, ^] and it corresponds to the 
upper branch of the Fundamental Diagram on Fig. 4. 

(b) For any moment of time the are infinitely many particles having gaps smaller than v (and 
hence zero normalized local velocities). Therefore at least for spatially periodic configurations 
we can apply Lemma 5.2 which guarantees that only gaps smaller than 2v will survive with 
time. Thus to study asymptotic properties it is enough to consider configurations having only 
two types of gaps: smaller than v and between v and 2v. 

Denote by X(L, m, n) the subset of admissible configurations x G X being spatially periodic 
with the spatial period of length L G R + , which contains exactly m G Z + particles with gaps 
belonging to the interval [0, v) and n G Z + particles with gaps belonging to the interval [v, 2v). 
Obviously p(x) = (m + n)/L. The set X(L,m,n) is invariant under dynamics (each time when 
the size of a gap crosses the threshold v one "small" gap becomes large and one "large" gap 
becomes "small") which immediately yields the exact value of the average velocity V(x) = -^h- 
On the other hand, by definition mv + n2v > L since the corresponding gaps fill in the segment 
of length / and lengths of both types of gaps are smaller than v and 2v respectively. Therefore 
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(p(x)L + n)v > L and hence n > L/v — p(x)L, which gives the lower bound 



nv nv L/v-p(x)L 

V(x) = = . — > v — — -. \ \ = I pix) - v. 

v ' m + n p(x)L p(x)L " y ' 

Observe, that choosing "small" and "large" gaps of length v — e and 2v — e for < e < 1 we 
see that the lower bound can be "almost" achieved. 

The upper bound of the average velocity in the hysteresis phase (i.e. when ^ < p(x) < i) 
follows from the existence of configurations with equal gaps of size larger than v for all densities 
from this segment. For the case p(x) > l/v the upper bound is calculated using the opposite 
length estimate nv < L. Then we get 

„, , nv nv L 

V(x) = = ——- < — — - = 1 p(x), 

v ' m + n p(x)L p(x)L 

which agrees with the weak normalization case. 

It remains to show that the region H is filled in densely by the pairs (p, V) corresponding to 
admissible configurations. To this end one considers all possible choices of the integer param- 
eters n,m and lengths of the corresponding gaps to get the result. Indeed, Vp G (^-, -) there 
exists an arbitrary large L such that pL G Z+. Choosing now various available combinations 
of positive integers m, n for which m + n = pL we can approximate V with the accuracy 

m + n pL 

□ 

Remark. By Theorem 5.1 for a.e. pair (V, p) G H there exists an admissible configuration 
x G X such that p(x) = p and V(x) = V. On the other hand, it might be possible that for 
some configurations having densities belonging to the hysteresis region the average velocity is 
not well defined and we claim only that all limit points of finite time velocities belong to the 
vertical segment corresponding to the given density. 

Corollary 5.3 Let x{r) G X(r), r > and p{x{r)) be well defined and let \/i,t v\ = v. Then 
V(x) = v if p(x(r)) < 2v +2r an d otherwise for a.e. point (p, V) in the curvilinear region 

H := {(p, V) : max(-^- - v, 0) < V < min(-^-, v)} 

there exists a configuration x(r) G X(r) with p(x(r)) = p,V(x(r)) = V, i.e. the region H 
corresponds to the hysteresis. 



6 Local velocities of both signs 

A close look to the previous analysis shows that we practically did not use the property that 
all particles move in the same direction, i.e. that P{v\ > 0) = 1. Now we explain the changes 
necessary to study this more general case. Consider an infinite configuration x(r) G X(r) and 
again interpret the values {vj}^ (which now may have both positive and negative signs, but 
still assuming that \vf\ < v) as local velocities for particles in the configuration x l (r). 

The presence of particles moving in opposite directions leads to a serious modification of 
the inequalities describing the violation of the admissibility condition for the i-th local velocity. 



20 



Actually this is the main and the most serious change comparing to the case of nonnegative 
velocities. Now we need to take into account not only the position of the succeeding particle, 
but also its velocity, as well as the corresponding quantities related to the preceding particle. 
In this more general case the i-th local velocity does not break the admissibility condition if 
and only if 

max(x*_ 1 (r), x*_ 1 (r) + v *_ : ) + r < min(x*(r), x\{r) + vf) — r 

< max(x*(r),x*(r) + vf) + r < min(x- +1 (r), x\ +1 (r) + v\ +1 ) — r. 

If for some i G Z and j G {i — 1, % + 1} the corresponding inequality is not satisfied we say that 
there is a conflict between the i-th particle and the j-th one and one needs to resolve it. In 
terms of gaps A* between particles the inequalities above can be rewritten as follows: 

A* > max(^, - v) +1 , v) - v) +1 ), j £ {i — 1, i} (6.1) 

Since the dynamics again will depend only on the sequence of gaps {A*} between particles, for 
each r > one can make the invertible change of variables (1.1) (described in the Introduction) 
to the case of 'point' particles with r = which we shall study further. 

Exactly as in Section 1 the strong normalization means that we reject (nullify) all velocities 
leading to a conflict, i.e 

Kf/t t\ _ j v t if (6-1) holds 
jV,l " i,J;J :_ \0 otherwise. 

The situation with the weak normalization is more delicate. The way how it was defined in 
Section 1 can be characterized as the only non-anticipating procedure allowing conflicting parti- 
cles to move simultaneously whenever possible. Following this idea we say that a normalization 
is weak if the positions of particles at the next time step rr- +1 := x\ + N w {v\, x l ) satisfy the 
conditions: 

xt+le ({xj + vf} if (6.1) holds 
* \{x f j, x 1 ^ 1 } if 3 a conflict of the particle i with the particle j — i ± 1. 

The 1st line describes the case when the admissibility condition holds, while the 2nd line shows 
what happens if it breaks down. Namely, if the i-th particle moves in the same direction as the 
j-th one then (by the non-anticipation property) the former assumes the previous position of 
the latter (x' +1 = #*■), otherwise the positions of the conflicting particles at time t + 1 coincide. 
The latter fact is the most important property here. 

If directions of all instant local velocities coincide then (6.2) defines the normalization 
uniquely. However if their signs are different then (6.2) implies only that 

x l + l = x] +l e [x% x]] n [x\ + v\, x) + uj]. 

Thus the set of weak normalizations is quite broad, for example it includes a random normal- 
ization when two mutually conflicting particles moving in opposite directions meet at a random 
point belonging to the segments described above. One can give a "natural" specific construction 
of M w normalizing local velocities in such a way that positions of particles at the next moment 
of time will be the same as if the particles would move simultaneously at continuous time with 
the given local velocities until the admissibility condition breaks down: 



v\ if (6.1) holds 

-A*_! if A*_i < -vj, v\ < 0, v\_ x < 

A^ if A\<vl vj>0, v*>0 



A 



ii u i ^ u i+l 

x v\ if A*_i < v\_ x - vf v\ < 0, v\_ x > 
^J- x v\ if A* < v\ - vl +1 , v\ > 0, v\ +l < 0. 
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After this long discussion of the definition of the normalization procedure it is surprising to 
find that all arguments used in the analysis of the case of positive velocities remain valid with 
only very slight changes. 

Lemma 6.1 The upper/lower densities p±(x t ) are preserved under dynamics. 

Proof. One uses the same estimates as in the proof of Lemma 2.3 except that now 2 particles 
may simultaneously leave or enter a given spatial segment / (instead of 1). Thus the total 
change of the number of particles in I is less or equal to 2 and hence 

\p(x\l)-p(x t+ \l)\-\I\<2. 

□ 

Lemma 6.2 Let x e X then \V(x,j, t) — V(x, i, t) \ a.s. Vi, j E Z. 

Proof. Again one follows the same argument as in the case of nonnegative local velocities. 
The only difference is that in the analysis of the connection between A* and A* now one needs 
to consider new cases related to negative local velocities. 

Additionally here instead of the uniquely defined weak normalization we need to consider 
an arbitrary one. If both v\ and vj +1 are nonnegative we are in the situation considered in 
Section 2. Therefore the cases (a) and (b) hold automatically. Nevertheless we formulate all of 
them to prove that A* > A* Wt G Z : 

(a) the condition (6.1) holds. Then obviously the argument used in Section 2 woks. 

(b) v\ > A', vj +1 > 0. Again one uses the same argument as in Section 2. 

(c) v\ < -A*_ r Then N w (y\,&) < M w {y\,^) < and Af w (vl +1 , x*) > M w { i>* +1 ,£*). Hence 

A< +1 = A* - M w {v\, x») + AfM + i, &) 

> A* - N w (v\, x*) +Af w (vl +1 , x*) = A? 1 . 

(d) vl>0, v\ +1 < and v\ - v* i+1 > A*. Then by definition A^ +1 > = A\ +1 . 

In the strong normalization setting one also considers the same cases and proves by induction 
that A* +1 > A* +1 - 2v (instead of . . . — v in the situation v\ > 0). New cases are the following 

(c') v\ < -A*_ x . Then 

N s {v\, X*) = v\ < -A*_! = M s (v\, X 1 ) 

and 

ACK* +1 ,x*)-A4(^ +1 ,^)>-2^ 
by the induction assumption. Hence 

A* +1 =%-K(vls*)+K(vl 1 ,x t ) 

>A\-2v- N s {v\, x l ) - Afsivj^x 1 ) + 2v 
= A* +1 . 
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(d') v\ > 0, v\ +l < and A* < vj - vj +l < A*. Then 

A* +1 = A* + w « - > A* + A* 
> -2-u + A* = -2w + A^ +1 . 

(d") u| > 0, < and A* < v\ - v\ +1 > A*. Then 

A^ +1 = A • > A • - 2-u = A* +1 - 2u. 

Note that the difference A^ +1 - A* +1 = -2v may be achieved only in the case (d'). The 
continuation of the proof is exactly the same as in Section 2, except for the change of 2v to Av 
in the last inequality. □ 

Using these results and applying exactly the same arguments as in the proof of Theorem 4.1 
one gets the uniqueness of the average velocity. 

Theorem 6.1 In the weak normalization case the set of limit points as t — > oo of the sequence 
{V(x,t)} te z depends only on the density p{x). 

t-i 

Theorem 6.2 (Fundamental Diagram) In the deterministic setting V (x) = lim ) ^ min(l/p, Vq). 

s =o 

Proof. Since at each moment of time t G Zo the local velocities of particles coincide, the 
condition (6.2) implies that 

x l +1 £ { x i + ^0) X \±l\- 

Thus the construction used in the proof of Theorem 4.2 remains valid in this case as well. □ 



7 Generalizations and Discussion 

7.1 Anticipating normalization 

Throughout the paper we consider only non-anticipating normalizations. In principle one might 
try to consider an anticipating normalization allowing at time t the i-th particle to move up 
to the position of the (i + l)-th particle rr*+} at time t + 1 rather than to x\ +1 . From the 
first sight this makes the normalization scheme more flexible. Unfortunately the anticipating 
normalization is not well posed since it turns out to be nonlocal. Namely a single change in 
the sequence of local velocities (say of the i-th one) may drastically alter the behavior of the 
system for particles having indices arbitrary far from the changed one (i.e. for j <C i). 

7.2 One-sided particle densities 

The density of a configuration in the way how it was defined in Section 2 depends sensitively 
on the statistics of both left and right tails of the configuration. A close look shows that in fact 
if all particles move in the same direction, say right, one needs only the information about the 
corresponding (right) tail, which allows to expand significantly the set of configurations having 
densities and for which our results can be applied. 

For a configuration x G X by a one-sided particle density we mean the limit 

p(x) :=limp(x,[0,£}). (7.1) 
The upper an lower one-sided densities correspond to the upper and lower limits. 
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Theorem 7.1 Let v\ > V«,t. Then all results of Lemma 2.1 and Theorems 4-1, 4-®> 5.1 
remain valid if one replaces the usual particle density p to the one-sided density p. 

Proof. The key observation here is that the assumption v\ > Vi, t implies that the movement 
of a given particle in a configuration x l G X depends only on particles with larger indices. 
Therefore if one changes positions of all particles with negative indices the particles with positive 
indices will still have the same average velocity. On the other hand, by Lemma 2.5 the average 
velocity does not depend on the particle index. This allows to apply the following trick. 

For each configuration x G X of density p{x) we associate a new configuration x G X defined 
by the relation: 



Then obviously p(x) = p(x) = p(x). 

Therefore for all purposes related to the average velocities all results valid for the configu- 



Note however that this trick does not work for the case of local velocities of both signs (con- 
sidered in Section 6), nor in the passive tracer analysis (Section 7.4). In both these situations 
statistics of particles with negative indices cannot be neglected. 

7.3 Nagel-Schreckenberg traffic fiow model 

The celebrated Nagel-Schreckenberg traffic flow model introduced in [16] for the lattice case is 
very similar to our case but additionally to the lattice setting it uses a bit different dynamics. 
In our terms this model differs from the main model introduced in Section 1 by that at each 
time step the previous normalized local velocity of the i-th particle is increasing by < a < a\ 
until it reaches v. One can think about a\ as an acceleration under the action of a (random) 
force (see e.g. [6]). Nevertheless the formalism elaborated in the present paper allows to study 
the continuum version of the Nagel-Schreckenberg model as well. In particular, in the weak 
normalization case one applies basically the same arguments as in Sections 2,3 and 4 since 
the distance between pair members cannot exceed C(v,a) < v 2 /a. Note however that the 
average velocity should be calculated in a more complicated way. Observe also that one can 
consider random accelerations of both signs a\ G (— oo, —a] U [a, oo) which makes the model 
more applicable. 

Mathematical formalism developed in the present paper can be applied with minimal changes 
to a number of other traffic flow models (discussed in detail, e.g. in a recent review [15]) allow- 
ing not only to study their continuum versions but also to get rigorous results in the original 
lattice setting which are absent at present. 

7.4 Passive tracer 

Following the idea introduced in [5] we study the dynamics of a passive tracer in the flow of 
particles imitating a motion of a fast pedestrian in a slowly moving crowd of people. 

Consider a pure deterministic setting {y\ = v) with the weak normalization and let T^x 
describe the flow of particles. The passive tracer occupies the position y l G H. at time t and 
moves all the time in the same direction. Before carrying out the next time step of the model 
describing the flow of particles, the tracer moves in its chosen direction to the closest (in this 
direction) position of a particle of the configuration T l v x. After that the next iteration of the 
flow occurs, the tracer moves to its new position, etc. 




ration x remain valid for x as well. 



□ 
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To be precise, let us fix a configuration x G X with p(x) > and introduce the maps 
: H. — > R defined as follows: 



r+y := min{xi : x t > y}, r x y := max{xi : x { < y}. 

Then the simultaneous dynamics of the configuration of particles (describing the flow) and the 
tracer is defined by the skew product of two maps - the map T v and one of the maps r. 1 * 1 , i.e. 

(x,y) ->T±(x,y) := (T v x,T±y), 

acting on the extended phase space IxR. The sign + or — here corresponds to the motion 
along or against the flow. We define the average (in time) velocity of the tracer 

V tl (x,t) :=(y'-y )/*, 

i.e. the total distance covered by the tracer (which starts at position y° G R) up to time t G Z + 
with the positive sign if the tracer moves forward, and the negative sign otherwise. 

Theorem 7.2 Let v\ = v Vi,t, Af = Af w , x G X and let x° i+l > xf Vi G Z. // the tracer 
moves along the flow (i.e. in the case T + ) then 

t^o? tt ( x _ fv if < p(x) <l/v 



V tI (x,t)^V(x) 



\l/p(x) othewise 



If the tracer moves against the flow (case 71 ) then V tI (x,t) V(x) — l/p(x). 

Proof. The assumption x Q i+1 > x® Mi G Z implies that x\ +1 > x\ Vi, t which allows to avoid a 
pathology related to the presence of several particles at the same position. In such a situation 
the tracer may "jump" through all of them in one time step. This cannot happen if r > in 
distinction to the case of point particles (r = 0). 

In the case of T + the tracer will run down one of the particles in the flow and will follow it, 
but cannot outstrip. Thus V tr (x,t) < -^5 V(x). 

Consider now the case when the tracer moves backward with respect to the flow, i.e. 71. 
Each time when the tracer encounters a particle, on the next time step this particle moves in 
the opposite direction and never will interfere with the movement of the tracer. Thus during 
time t > the tracer meets exactly t particles which gives 

(-V tr (x,t) + V(x,t))tp(x) = t. 

Therefore 

V tI (x,t) = -l/p(x) + V(x,t). 

□ 



Using similar arguments in the case of the strong normalization one can show that V tT (x, t) in 
the gaseous phase of the particle flow has the same asymptotic as in the weak normalization case. 
Since the flow in the fluid phase demonstrates hysteresis the same phenomenon is unavoidable 
for the passive tracer as well. 
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7.5 Multidimensional generalization 

The constructions used in this paper are essentially one-dimensional. Still at least some direct 
generalizations are possible. Let x\ G M d , d G Z + and denote by (xf)j the j-ih coordinate of 
the ci-dimensional vector x\. We say that a configuration x*(r) is admissible if 

max((x-(r))o) + r < min((x* +1 (r)),) - r Vi G Z. (7.2) 

All results of Sections 2, 3, 4.1, and 6 hold in this setting. Unfortunately the assump- 
tion (7.2) implies that a natural multidimensional generalization of the notion of density of 
the configuration x l (r) turns out to be equal to zero for any admissible configuration. How- 
ever densities for one-dimensional projections are well defined and for them the Fundamental 
Diagram type results are readily available. 

7.6 Open problems and conjectures 

Our construction give a very precise information about the asymptotic properties of DDS under 
consideration in the deterministic setting. In the random setting we prove only the uniqueness of 
the average velocity. From the results of Section 2 it follows that the mathematical expectation 
of lower/upper average velocities are well defined but we are not able to calculate them. On 
the other hand, we can formulate a conjecture that the limits as time goes to infinity of finite 
time average velocities are deterministic. In other words, the Law of Large Numbers is valid 
for the sequence of finite time average velocities. 

An important question is whether the dynamical coupling of pairs of processes with equal 
densities under the weak normalization is nearly successful. Let V be the common distribution 
of the i.i.d. local velocities. As we know in the pure deterministic setting when the distribution 
V is concentrated at a single point {v} the dynamical coupling needs not to be successful. 
Nevertheless we conjecture that for each nontrivial distribution V the nearly successful coupling 
takes place. Moreover, the non-triviality of the distribution V should lead to the existence and 
uniqueness of the translationally invariant measure of the Markov chain described by the DDS. 
Proofs of results of this sort need the development of an additional probabilistic apparatus and 
will be discussed elsewhere. 
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